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We reveal the significance of kinetic-driven multiple-spin interactions hidden in geometrically- 
frustrated Kondo lattice models. Carefully examining the perturbation in terms of the spin-charge 
coupling up to the fourth order, we find that a positive biquadratic interaction is critically enhanced 
and plays a crucial role on stabilizing a spin scalar chiral order near 1/4 filling in a triangular lattice 
case. This is a generalized Kohn anomaly, appearing only when the second-order perturbation is 
inefficient because of the degeneracy under frustration. The mechanism is potentially common to 
frustrated spin-charge coupled systems, leading to emergence of unusual magnetic orders. 



Localized spins in metal interact with each other via 
effective interactions mediated by the kinetic motion of 
itinerant electrons. An underlying mechanism is theoret- 
ically established as the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interaction The precise form of the inter- 

action sensitively depends on the electronic structure of 
the system as well as the coupling between charge and 
spin. This offers the possibility of various magnetic or- 
derings in the spin-charge coupled systems. 

Recently, unusual magnetic orderings have been of in- 
tensive research interest in geometrically-frustrated spin- 
charge coupled systems. In particular, a spin scalar chi- 
ral ordering has attracted considerable attention, since 
it underlies numerous fascinating phenomena, such as 
the unconventional anomalous Hall effect [3, El- As a 
prominent example, a scalar chiral ordering with a four- 
sublattice noncoplanar spin configuration was discussed 
for the Kondo lattice model on a triangular lattice [B-Q ■ 
The particular order appears in two regions, near 3/4 fill- 
ing and 1/4 filling 0. While the former was predicted on 
the basis of the perfect nesting of the Fermi surface @ , 
the latter is unexpected from the nesting scenario. Sur- 
prisingly, the 1/4 filling phase is stable up to much larger 
spin-charge coupling than the 3/4 filling one Q ■ The ro- 
bustness of the noncoplanar order near 1/4 filling casts a 
fundamental question on the mechanism of chiral order- 
ing. 

In this Letter, we point out the significance of kinetic- 
driven multiple-spin interactions for the noncoplanar or- 
dering under geometrical frustration. Carefully exam- 
ining the perturbation in terms of the spin-charge cou- 
pling up to the fourth order, we find that the second- 
order RKKY contributions are degenerate among several 
different magnetic orderings because of the frustration, 
whereas the degeneracy is lifted by fourth-order contri- 
butions. Among several terms, a positive biquadratic 
interaction is critically enhanced and plays a crucial role 
in the emergence of the noncoplanar order. We ascribe 
the mechanism of the critical enhancement to a Fermi 
surface effect a generalized Kohn anomaly at a higher 
level of perturbation. It opens a local energy gap at sev- 



eral points on the Fermi surface, which develops continu- 
ously to a full gap as increasing the spin-charge coupling 
and approaching 1/4 filling. This naturally explains the 
robustness of the chiral ordered phase. 

We consider the Kondo lattice model on a triangular 
lattice, following the previous studies @, 0, 0, E3] • The 
essential result is common to other frustrated lattices, as 
mentioned later. The Hamiltonian is given by 



H = -t > ' (c| a Cj >a + h.c.) - J n y ] cT a a a pCi^ ■ S, : 



(i) 



where c\ a (cj )ai ) is a creation (annihilation) operator of 
conduction electron at site i with spin a, <x is the Pauli 
matrix, and Si is a localized moment. The sum in the 
first term is taken over the nearest-neighbor sites on the 
triangular lattice. The sign of the spin-charge coupling 
Jh does not matter, since we assume Si to be a classical 
spin. 

Given a specific spin configuration {S^}, it is possi- 
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FIG. 1. (Color online) Ordering patterns of localized spins 
considered in the perturbation calculations: (la) ferro- 
magnetic, (2a) two-sublattice collinear stripe, (3a) three- 
sublattice 120° noncollinear, (4a) four-sublattice all-out, and 
(4f) four-sublattice coplanar 90° flux orders. The right panel 
is the phase diagram of Hamiltonian fTJ near 1/4 filling. 
Arrows (i) and (ii) indicate the parameter changes used in 
Fig. Hfd). PS represents a phase-separated region and the 
bold dotted line at 1/4 filling denotes an insulating (4a) state. 
See also Figs. 1 and 3 in Ref. 0- 
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ble to obtain the exact free energy F = — 7 1 ^ a log(l + 
cxp( — (3E a )) from the eigenvalues of H, {E a } (T is tem- 
perature and [3 = 1/T). In the previous study, from 
the comparison of the ground state energies, two of the 
authors discovered a scalar chiral phase with a four- 
sublattice noncoplanar spin order [Fig. Qj4a)] near 1/4 
filling Q ■ As shown in the right panel of Fig. [1] the chi- 
ral order sets in by an infinitesimal Jh from slightly less 
than 1/4 filling, and extended to a large Jh insulating 
region at 1/4 filling Q . 

In order to elucidate the origin of the robust noncopla- 
nar ordering, we consider the perturbation in terms of 
Jh- The analysis provides coherent understanding of the 
phase diagram not only in the small Jh region but in the 
large Jh insulating region, as we will discuss later. Tak- 
ing the second term in eq. ([l} as the perturbation {%') to 
the first term (%o)i we use the standard perturbation the- 
ory to calculate the energy for various ordering patterns. 
We here consider a collection of ordered states up to four- 
sublattice ordering, which appeared in the phase diagram 
in the small Jh region; namely, (la) ferromagnetic, (2a) 
two-sublattice stripe, (3a) three-sublattice 120° coplanar, 
(4a) four-sublattice all-out, and (4f) four-sublattice flux 
orders, as shown in Fig. [T] The free energy is expanded 
in the form 



F — Fq = -Tlog(Vexp(- J U'(T)dT 
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where Fq is the free energy obtained from %q 1 t is the 
imaginary time, and T represents time-ordered product. 
(• • • }con means the average over the connected diagrams. 
From the spin rotational symmetry, odd order terms do 
not appear in the expansion. 

By decoupling the brackets (• ■ • ) con with the Wick's 
theorem, each term in eq. @ can be written in a 
compact form. The second-order term gives the well- 
known RKKY form [1^; F^ = - J 2 £ q |S q | 2 X o(q), 
where the bare magnetic susceptibility is defined as 
Xo(q) = -iEk{/f( £ k-q) - /F(£k)}/(ek- q - £k), and 
S q = jjJ2jSje t<i rj (If is the Fermi distribution func- 
tion, ek is an eigenvalue of T~lo in the momentum basis, 
and N is the number of sites). Specifically, for each mag- 
netic order, J 1 *- 2 ' can be rewritten as 

f X o(0,0) for (la) 

F& = —J 2 x < Xo (27r/3,27r/3) for (3a) (3) 
[xoKtt) for (2a) (4a) (4f). 

Here, we take the triangular lattice in the form of a square 
lattice, as shown in inset of Fig. [^a). Note that the (2a), 



(4a), and (4f) states have the same energy within the 
second-order perturbation. 

The fourth-order term in eq. ([2]) is calculated in a 
similar manner. For example, the result for the four- 
sublattice orders is summarized into a form of effective 
multiple-spin interactions as 

^ (4) = |(t) 4 £E E E 

v 7 cj p k ji j2,i3,i4 = l~4Q,.Q m ,Q„ 
e iS Jl -(Qi+Q m +Q n ) e -iSj 2 -Qi e -iS j3 ^Q me -iS j ^Q„ x 

G k (o;p)G'k + q ?i ( w p)G k+ Q m+ q n ( UJ p) < ^k+Q l +Q m +Q n ( w p) 
x [(S h ■ S, 2 )(S J4 • S js ) + (S h x S j2 ) • (S j4 x Si,)] 

= j h J2 [ JSi • s j + B ( Si ■ s ^ 2 + B *( s * x s j) 2 ] 

l<i<j<4 

+ [ M ( s -i ■ s j)( s * ' s k) + M *& x Si) ■ (S, x S fc )] 
+ [l(SrS J )(S t -S ( )H*(S,xS J ).(S fc xS ! )]. 

(4) 

Here the sum of Q; i?Tt , n runs over (0,0), (7r,0)(= Q a ), 
(0, 7r)(= Q b ), and (tt,tt)(= Q c ) (Qa, Q&, and Q c arc 
four-sublattice ordering wave vectors [6[ ) . The indices i, j 
denote the four sublattices, and 5j = (0,0), (1,0), (0,1), 
and (1, 1), respectively [see inset of Fig. [2fa)]. G k (wp) = 
{iu>p — (e k — A 4 )} i s the noninteracting Green's function, 
where lu p is the Matsubara frequency and fi is a chemical 
potential. Thus, the fourth-order contribution includes 
three-spin and four-spin effective interactions in addition 
to two-spin interactions. Each coefficient is given by the 
sum of four Green's functions products. 

First let us discuss the result of second-order pertur- 
bation. The second-order term is proportional to \o, a s 
shown in eq. ([3]). The comparison of xo as a function 
of the electron density n is presented in Fig. (Ha) . Since 
the magnetic order with largest xo is stabilized when Jh 
is turned on, the result indicates that the ferromagnetic 
order appears for 0.0 < n < 0.15 and 0.8 < n < 1.0, and 
the three-sublatticc 120° coplanar order is realized for 
0-3 < n < 0.55. See also Fig. 3 in Ref. 0. Meanwhile, 
Xo(7r,7r) is the largest in the remaining regions, in which 
we have to consider the next fourth-order terms to lift 
the degeneracy between (2a), (4a), and (4f) states 
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The fourth-order free energies in the remaining re- 
gions for (2a), (4a), and (4f) are indicated in Figs. (5Jb) 
and [UJc). The results show that (2a) two-sublattice 
stripe phase has the lowest fourth-order contribution in 
a wide range of n, while (4a) four-sublattice all-out or- 
der becomes lowest in narrow windows at n ~ 0.225 and 
n ~ 0.75. The stability of (4a) state at n ~ 0.75 is ex- 
plained by the perfect nesting of the Fermi surface at 
3/4 filling @. On the other hand, the reason for the 
stabilization at n ~ 0.225 is not obvious. 



3 



( a ) 0.6 
0.5 
0.4 
%D 0.3 
0.2 
0.1 



0.0 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

n 



(b)o.30 






0.25 


(2a) 


1 


0.20 






^0.15 
I3O.IO 


(4f) 




sT~ 0.05 
0.00 




A 






-0.05 






-0.10 






-0.15 





0.18 0.20 0.22 0.24 0.26 0.28 

n 



Ho 

8 




(2a) j \ 


J— » 


(4a) ,\ 

< 4f ) A 


£~~ 2 


-2 






-4 





0.64 0.68 0.72 

n 



076 



FIG. 2. (Color online) (a) Xo(Q) i n eq. (O as a function of 
n. The inset represents a square lattice with diagonal bonds 
which is topologically equivalent to the triangular lattice, (b), 
(c) The fourth-order free energies of (2a), (4a), and (4f) states 
in the region near (b) n ~ 1/4 and (c) n ~ 3/4. In the 
calculation, we take 800 x 800 grid points in the extended 
first Brillouin zone and T = 0.03t. 
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FIG. 3. (Color online) Coefficients of the different terms in 
the fourth-order effective Hamiltonian [eq. ©] as a function 
of the electron density n. 



In order to gain an insight into how the (4a) state is 
stabilized at n ~ 0.225, we carefully analyze the fourth- 
order contributions in cq. ^ term by term. We plot 



the coefficient of each term, the product of the Green's 
functions, as a function of n in Fig. |3l The results 
show that all the coefficients show drastic changes around 
n ~ 0.225: In particular, the coefficient of biquadratic in- 
teraction, B, shows a pronounced peak with the largest 
contribution. This indicates that the effective positive bi- 
quadratic interaction plays a dominant role at n ~ 0.225. 
In general, the positive biquadratic interaction favors a 
noncollinear spin configuration. Particularly, among the 
degenerate (2a), (4a), and (4f) states, the (4a) noncopla- 
nar all-out order is naturally selected by this interaction. 

The emergence of the positive biquadratic interaction 
is remarkable, because the effective biquadratic interac- 
tion has a negative coefficient in most cases. For ex- 
ample, in the Hubbard-type Hamiltonian, a fourth-order 
perturbation from the strong-coupling limit usually leads 
to a negative biquadratic interaction. Moreover, in lo- 
calized spin systems, a spin-lattice coupling and ther- 
mal/quantum fluctuations also give rise to negative bi- 
quadratic interactions. Our result indicates that an un- 
usual positive biquadratic interaction is induced by the 
kinetic motion of electrons at a special filling. 

Why is the positive biquadratic interaction enhanced 
here? In order to elucidate the origin, let us look at the 
electronic state at n ~ 0.225 in detail. FigureHJa) shows 
the Fermi surface at n = 0.225, which is circular, with no 
substantial nesting tendency. At first sight, this feature- 
less Fermi surface does not lead to any type of singularity. 
It, however, conceals an unexpected instability toward 
four-sublattice ordering. The key observation is that the 
ordering vectors of four sublattice order, Q a = (jr, 0), 
Qb = (0,7r), and Q c = (tt,tt), connect the particular 
points on the Fermi surface when n reaches ~ 0.225; the 
tangents of the Fermi surface are parallel between the 
connected pairs — see Fig. Ufa). These connections give 
rise to a singularity analogous to the 2k f -singularity or 
Kohn anomaly, observed in an isotropic electron gas 12 1 . 
It is worthy to note that, in contrast to the Kohn ano maly 
usually discussed in the second order perturbation [12} . 
the singularit y co nsidered here arises in the fourth-order 
perturbation [13} . 

The scenario is confirmed by numerically calculating 
the contributions only from the parts of Fermi surface 
connected by Q a , Qb, and Q c . We calculate the coef- 
ficients of the effective four-spin interactions in cq. @ 
by limiting the k summation within the small hatched 
patches shown in Fig. HJa) . As shown in Fig. HJb) , the 
contributions from the small patches with 6 = ir/40 lead 
to a positive and largest biquadratic interaction B, which 
well accounts for the result in Fig. [3J This indicates that 
the anomalous behaviors in Fig. [3] can be attributed to 
the particle-hole process within such small area in the 
momentum space. 

The critical enhancement is, in fact, a divergence 
in the limit of T — > in the perturbative calcula- 
tions. The T dependence is analytically calculated by 
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FIG. 4. (Color online) (a) The Fermi surface at n = 0.225 
in the extended Brillouin zone scheme. The connections 
of the Fermi surfaces by the four-sublattice wave vectors 
Qa = (tt, 0), Qt = (0, tt), Q c = (tt, tt) are shown. Six hatched 
squares in the first Brillouin zone indicate the patches used in 
the calculations of the data in (b). The patches are 5 xS each, 
and centered at the special k points, (±7r/2,0), (0,±7r/2), 
and (±7r/2, =p7r/2), which are connected with each other by 
the four-sublattice wave vectors, (b) Coefficients of the effec- 
tive four-spin interactions in eq. ((4]) obtained by taking the 
k summation only within the patches in (a) with 5 — tt/40. 

(c) The peak values of B = T £ k £„ p [G£(o, p )] 2 [G° k+Q (oj p )] 2 
with patch S = tt/S as a function of T. The peak value of 
B calculated by the integration over the entire first Brillouin 
zone is also plotted. Two straight lines are the fitting to T~ 1,5 . 

(d) The Fermi surface evolution in the (4a) phase along the 
arrows in the right panel of Fig. [1] (i) from Jh = to 0.75 
(every 0.25) at a fixed n = 0.2255 and (ii) from n = 0.231 
to 0.249 (every 0.006) at a fixed Jh = 0.75. Only the first 
quadrant of the folded Brillouin zone is shown. 



approximating the connected fragments of the Fermi sur- 
faces by simple biquadratic functions centered at the 
connected k points. The dominant contribution to 
the divergence comes from the terms proportional to 
TE k E Wp [G°H)] 2 [G£ +Q H)] 2 (= B) in eq. ©; it 
leads to quick divergence oc T -15 , reflecting the fact 
that the tangents of the Fermi surface are parallel to 
each other at the connected points. The details of the 
calculaions will be reported elsewhere. Indeed, the nu- 
merically calculated values of B as well as B increase 
with oc T~ 15 as T 0, as shown in Fig. HJc). 

The divergence leads to a local gap formation at the 
connected points on the Fermi surface. FigureHJd) shows 
the Fermi surface evolution in the (4a) phase, as we 
change the parameters (n, Jh) along the arrows (i) and 
(ii) on the right panel of Fig.[T] As soon as Jh is switched 



on, the Fermi surface is modified at the connected points 
reflecting the local gap opening. The gap extends to 
shrink the Fermi surface as Jh increases and n becomes 
closer to 1/4 filling; eventually, the Fermi surface dis- 
appears in the insulating state at 1/4 filling. This con- 
tinuous change clearly demonstrates that the generalized 
Kohn anomaly revealed by the perturbative argument 
provides an essential ingredient for the robust chiral or- 
dering including the insulating region. 

Let us emphasize the importance of the geometrical 
frustration in this generalized Kohn anomaly. Frustra- 
tion leaves degeneracy at the level of second-order pertur- 
bation, and prepares the stage for the fourth-order term 
to select the type of magnetic order. In this sense, ge- 
ometrical frustration is an indispensable prerequisite for 
the positive biquadratic interaction to play an active role. 
In fact, when the frustration is reduced by introducing a 
spatial anisotropy in the triangular lattice structure, non- 
coplanar order is not stabilized within the four-sublattice 
order, even though the positive biquadratic interaction 
is present in the fourth-order terms: the effect of posi- 
tive biqadratic interaction is completely masked by the 
second-order RKKY interaction. Meanwhile, for other 
frustrated lattices, such as checkerboard, face-centered- 
cubic, and pyrochlorc lattices, we found that similar four- 
sublattice order is commonly stabilized near 1/4 filling by 
the positive biqiuadratic interaction [l4j|. This suggests 
that the present mechanism is robust widely in frustrated 
spin-charge coupled systems. 

To summarize, we have carried out the perturbation 
expansion up to the fourth order in Jn/i and constructed 
the effective multiple-spin interactions for the ferromag- 
netic Kondo lattice model on a triangular lattice. Among 
the kinetic-driven four-spin interactions, we have found 
that a biquadratic interaction takes a large positive coeffi- 
cient and increases divergently as T — > near 1 /4 filling. 
This signals an instability toward a noncoplanar chiral 
phase near the 1/4 filling with a local gap formation. 
The origin of large positive biquadratic interaction is as- 
cribed to the Fermi surface connection by the ordering 
wave vectors of four sublattice order, which we call the 
generalized Kohn anomaly. Our result provides compre- 
hensive understanding of not only the emergence but also 
the robustness of the chiral ordering. The mechanism is 
universal in frustrated spin-charge coupled systems, pos- 
sibly even when electron correlation between conduction 
electrons is included unless the Fermi surface is destroyed 
by electron correlation. 
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